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ABSTRACT 



We report the detection of four new extrasolar planets in orbit around the moderately active stars HD 63765, HD 104067, HIP 70849, 
and HD 125595 with the HARPS Echelle spectrograph mounted on the ESO 3.6-m telescope at La Silla. The first planet, HD 63765 b, 
has a minimum mass of 0.64 Mj„p, a period of 358 days, and an eccentricity of 0.24. It orbits a G9 dwarf at 0.94 AU. The second 
planet, HD 104067 b, is a 3.6 Neptune-mass-planet with a 55.8-day-period. It orbits its parent K2 dwarf, in a circular orbit with a 
semi-major axis of a=0.26 AU. Radial velocity measurements present a a;500-day-oscillation that reveals significant magnetic cycles. 
The third planet is a 0.77 Neptune-mass-planet in circular orbit around the K4 dwarf, HD 12595, with a 9.67-day-period. Finally, 
HIP 7849 b is along-period (5< P <75 years) and massive planet of m. sin; a;3.5-15Mjup that orbits a late K7 dwarf. 
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1. Introduction 

The HARPS planet-search program has been going on for 6 
years (since 2003) at the 3.6-m ESO telescope located at La 
Silla Observatory, Chile ( |Pepe et al.|[2000{ [Mayor eraLl[2003 1. 
Both HD 63765 and HD 104067 were part of the high-precision 
subprogram of the HARPS Guaranted Time Observions (GTO), 
which aimed at detecting very low-mass planets in a sample of 
solar-type stars already screened for giant planets at a lower pre- 
cision with CORALIE Echelle spectrograph mounted on the 1 .2- 
m Swiss telescope on the same site ( Udry et al. 2000 ). HIP 70849 
and HD 125595, on the other hand, are part of the HARPS 
volume-limited sample program of the HARPS GTO. The sam- 
ple is composed of about 1200 stars within 57.5pc, among which 
850 G, K dwarfs that are known as single stars with low activ- 
ity and low v. sin /. The goal of this program was to monitor the 
850 stars at the 3 ms"' precision level to detect all the gas giants 



Send offprint requests to: Damien Segransan, e-mail: 
Damien. Segransan@unige . ch 

• The HARPS radial velocity measurements discussed in this paper 
are only available in electronic form at the CDS via anonymous ftp to 
cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi- 
[biri/qcat?J/A+A/ 

^ Based on observations made with the HARPS instrument on the 
ESO 3.6 m telescope at La Silla Observatory under the GTO programme 
ID 072.C-0488. 



present in this statistically robust sample. Observed distributions 
of planets' orbital elements and masses, as well as host star char- 
acteristics should then be compared to theoretical predictions to 
improve our knowledge of planet formation and evolution. 

So far, HARPS has allowed the detection (or has contributed 
to the detection) of more than 85 extrasolar planet candidates 
(the late st batch of which was presented in a series of parallel pa- 
pers e.g. Forveille et al.|20lf Levis et al.|2011 Mordasini et al 



201 1 ; Moutou et al.|2011||Lo Curto et al.|2010t|Naef et al.|20Tof 
among the a!460 known today. But more specifically, HARPS is 
unveiling the tip of the iceberg of the very low-mass planet pop- 
ulation (see ' Lovis et al.|2009 Mayor et al.|200 9') which when 
combined with forthcoming results allow extrapolating results of 
their mass distribution. |Lovis et al.| ( [2009l l show that «30% of the 
stars surveyed in the HARPS high precision program have plan- 
ets with masses below 30-earth masses on periods shorter than 
50 days. Furthemore, a!80% of those planets are in multplanetary 
systems. HARPS' high stability not only gives access to super- 
Earth's at short periods but is also sensitive to Neptune-mass- 
planets at a longer peiiod, (e.g., HD 69830 b, Lovis et al. 2006|. 
The discovery of those intermediate-mass planets (1-10 MNept) 
in the 0. 1 to 0.5 AU separation range is of interest for planetary 
formation scenarios since those planets are hardly formed ac- 
cording to Ida & Lin ( 2005 1 or at least, should be rare according 
to |Mordasini et al.| ( |2009"^ 
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Fig.l. Left : Chromospheric activity index periodogram for HD 63765, HD 104067, HD 125595, and HIP 70849 . Right : Ca II H 
emission region. The large re-emission at the bottom of the Ca II H absorption line at /I = 3933.66 A is an indicator of chromospheric 

activity. 



In this paper we report the detection of four planets orbit- 
ing moderately active stars. Two of them are Neptunian planets 
at relatively short periods orbiting HD 104067 and HD 125595. 
The other two are gas giants at longer periods in orbit around 
HD 63765 and HIP 70849. The paper is organized as follows. In 
the first section we briefly discuss the primary-star properties. 
Data analysis is presented in section |3.1| while radial-velocity 
measurements and orbital solutions are discussed in section [3] 
In section l4j we provide concluding remarks. 

2. Stellar characteristics 

The ffective temperatures, gravity, and metallicities were derived 
from the spectroscopic analysis of the HARPS targets as de- 
scribed in Sousa et al. ( 2008) ). We used the improved Hipparcos 
astrometric parallaxes rederived by van Leeuwen ( 2007 1 to de- 
termine the absolute V-band magnitude using the apparent visual 
magnitude from Hipparcos (ESA 1997). Metallicities, together 
with the effective temperatures and My were then used to esti- 
mate basic stellar parameters (ages, masses) bases on theoretical 
isochrones from Girardi et al. (2000) and a Bayesian estimation 
method described in |da Silva et al.| (2006 r] Individual spectra 
were also used to derive both the bisector inverse slope ( BIS 



) of the HARPS cross-correlation function ( CCF ), as defined 



by Queloz et al. (2000), as well as a measurement of the chro- 



' The web interface for the Bayesian estimation of stellar param- 
eters, called PARAM 1.0, can be found at |http://stev.oapd.inaf.it/cgI^ 



mospheric activity index log (Rfjj^), following a similar recipe 
as used by 'Santos et al. ( 2000[ l for CORALIE spectra. Using 
theN oyes et al. ( 1984| l empirical calibration of the rotational pe- 
riod vs. activity, we estimated rotational periods. We also de- 



|bin/param| 



rived the v. sin (;) from the HARPS spectra following the |Santos 
|et al. ( 2002[ i approach. The HIP 70849 spectral type is, however, 
too late (K7V) to carry out [Sousa et al.| ( [2008) 's analysis. For 
such a low-mass star, we have to rely on near-infrared photomet- 
ric measurements, as well as on accurate mass-luminosity rela- 
tions ( [Baraffe et al.|1998|[DeIfosse et al.|2000l l, to derive its fun- 
damental physical parameters. The Ust of all stellar parameters 
can be found in Table [T| 



The four stars have activity indexes, log [Ru^], ranging be- 
tween -4.69 and -4.80, and they show a large Ca II re-emission at 
A - 3933.66 A (see Fig. Ill), revealing significant chromospheric 
activity possibly induced by stellar spots or plagues. In addition, 
log [Rfj/A indexes display variability on different time-scales as 
observed on their periodograms (see Fig. fT). Short time-scale 
variation can often be matched to the rotation period of the star 
(such as for HD 63765 and HD 125595 ), while longer time 
scale variability in the chromospheric activity index could betray 
the existence of stellar magnetic cycles (such as for HD 63765, 
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HD 104067 and HD 125595 ). HIP 70849 does not display any 
significant variabilitiy in the log [Rf]^) index. 

For some stars, a linear correlation between the activity index 
and the radial velocities is observed (see Fig.|3]l, leading the way 
to possible stellar activity detrending around moderately active 
stars. However, this topic is beyond the scope of this paper so we 
decided, instead, to estimate the overall stellar radial velocity 
jitter. To do so, we measured the radial velocity dispersion on 
a sample of stars with similar log [RffjA and spectral type that 
show no planetary. For each star, a radial velocity jitter of a few 
m.s~' was derived (see Table [Iji, which was used as an external 
white noise in the planet search algorithm (see Sect. |3.1[ ). 

3. Radial-velocity measurements and orbital 
solutions 

3.1. A Bayesian genetic algorithm witli Markof cliains 

Since 2003, we have been developing and using a genetic al- 
gorithm to retrieve multiple Keplerian orbital solutions (see 
Segransan] |2011[ in prep). The main advantage of genetic al- 
gorithm compared to other advanced methods such as "Markov 
chain Monte Carlo" (hereafter MCMC) with parallel tempering 
(for illustration see 'Gregory"2007 ) is that it allows probing the 
full model parameters space in a extremely efficient way (so- 
lutions of complex planetary systems solutions are derived in a 
few minutes instead of hours to run MCMC from scratch on a 
MacBook Pro). However, the population at the end of the evo- 
lution is not statistically reliable owing the intrinsic nature of 
genetic algorithm based of genome crossover and mutation. We 
therefore added a MCMC module with Metropolis hasting to 
overcome this problem. The algorithm is based on the formalism 
described in Colli er Cameron et al.| ( |2007] l with several thou- 
sands chains drawn from the final genetic algorithm's popula- 
tion. Each chain runs several thousand times to retrieve a sta- 
tistically reliable population. For poorly constrained planetary 
systems, the Bayesian formalism described in Gregory ( 2007| l 
is used. A full overview of the algorithm will be described in a 
forthcoming paper (see Segransan|2011| in prep). 



variation cr=3.41 m.s"\ yielding a reduced ;t'^ of 1.30. Figure |3] 
shows the HARPS radial velocities and the corresponding best- 
fit Keplerian model. The resulting orbital parameters for the 
first planet are P=358days, e=0.24 and /r=20.9ms"' , imply- 
ing a minimum mass m. sin (/)=0.64 Mjup and a semi-major axis 
fl=0.940 AU. Despite its proximity to one year, we are confident 
that the 358-day period is the correct one. Indeed, the genetic 
algorithm does not display any other significant solution in a 
period-eccentricity diagram (see Fig. l2]i. This is reasonably ex- 
plained by the large radial-velocity amplitude (K = 20.9 ms"') 
compared to the residual noise (cr=3.41 m.s"') and the fact that 
75% of the orbital phase is covered (see Fig. l3]l. In addition, 
the mcmc analysis shows that the probability distribution func- 
tions (hereafter PDF) of all parameters are gaussian with no 
obvious correlations between them as illustrated by the semi- 
amplitude vs. eccentricity diagram (see Fig.|2]i. Orbital elements 
for HD 63765 b are Hsted in Table [3] 



3.2. A one-year period Jovian planet around HD 63765 




16 17 18 19 20 21 22 23 
k1 [m/s] 

Fig. 2. Top:genetic algorithm population at the end of evolution 
represented in a period-eccentricity diagram. The dark curve is 
the maximum eccentricity a planet can reach without colliding 
the star. Bottom:joint probability density function of the radial 
velocity semi-amplitude and the eccentricity. The red, blue and 
cyan contour represent the 1, 2 and 3 cr confidence intervals of 
the joint PDF. 



HD 63765 was part of the CORALIE planet search survey (Udry 
et al. 20001 before being monitored at higher precision with 
HARPS. Since March 2000, 56 measurements have been ac- 
quired with CORALIE with a typical signal-to-noise ratio of 
46 (per pixel at 550 nm) and a mean measurement uncertainty 
(including photon noise and calibration errors) of 4.5 ms '. In 
June 2007, we upgraded the instrument that slightly changes 
its zero point (Segransan et al. 2010i We therefore adjusted 
an instrumental offset between the measurements taken before 
and after June 2007 (namely CORALIE-98 and CORALIE-07). 
Since December 2003, we have recorded a total of 52 mea- 
surements with HARPS with abetter signal-to-noise (SNRsil55) 

than CORALIE as well as a better final accuracy (60cm ', in- 3.3. A 55-day-period Neptune-like planet around HD 104067 
eluding photon noise and calibration errors). Despite the lower 
accuracy of the CORALIE spectrograph, we decided to include 
the CORALIE data because of their longer time span, mainly 
to increase the accuracy on the orbital period of the planet. 
Altogether, 108 radial velocities were gathered. A stellar veloc- 
ity jitter of 3.0 m.s"' was quadratically added to the mean radial 
velocity uncertainty. A single-planet Keplerian model was ad- 
justed to the data and produced residuals that show a level of 



HD 104067 has been observed with HARPS since February 
2004. Altogether, 88 radial-velocity measurements with a typical 
signal-to-noise ratio of ~ 172 (per pixel at 550 nm) and a mean 
measurement uncertainty (including photon noise and calibra- 
tion errors) of 43 cms"' were gathered. A stellar velocity jitter 
of 3.0 m.s"' was quadratically added to the mean radial velocity 
uncertainty. 
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Table 1. Observed and inferred stellar parameters for HD 63765,HD 104067, HIP 70849 and HD 125595. 



Parameters 



HD 63765 



HD 104067 



HD 125595 



HIP 70849 



Sp. T. 

V 

B-V 

J 

H 

K 

7T 



My 

[Fe///]'" 
log(g)<" 

Age 

Stellar litter 

V sin (0*^' 

Pro, 



[mas] 



G9V 

8.1 
0.75 



30.07 ± 0.56 



K2V 
7.93 
0.99 



47.47 ±0.90 



K4V 
9.03 
1.107 



35.77+ 1.04 



K7V 

10.36 

1.427 
7.639±0.023 
7.006+0.061 
6.790+0.027 
42.42+2.09 



[K] 
[dex] 

[Mo] 
[Gyr] 

[m/s] 

[kms-'] 
[days] 



5.49+0.04 
5432+19<" 
-0.16+0.01 

4.42+0.03 
0.865+0.029<-''> 

4.77+4.07^' 

-4.736+0.037*2) 

3.0 

1.63 
23 



6.31+0.05 
4969+72<" 
-0.06+0.05 

4.47+0.13 

0.791+0.020<^' 

4.33+3.99<3) 

-4.743+0.023<2' 

3.0 

1.61 
34.7 



6.80+0.06 

4908+ 87<" 

0.02+0.06 

4.32 +0.20 
0.756 + 0.017<'> 

3.53+ 3.59 '3' 

-4.787+0.037<2> 

2.5 

1.50 
40.1 



8.50+0.10 
4105 +130''*> 



0.63+0.03<'" 

I-5W 

-4.697<5' 

3.5 

1.93 



Notes. Near infrared photometry is taken from the 2M ASS All-Sky Catalog, jCutri et al.||2003[ |Skrutskie et al.| |2006"). As trometric data and 
visual photometry come from from Hipparc os Catalogs, l|ESA|1997 [van Leeuwen|2007|). Other parameters come from: (l) Sousa et al.H2008] l 
spectroscopic analysis, (2) HARPS CCF, (3) |Girardi et al.| ( |2000| ) models, (4) |Baraffe et al.| ( [T998l > models, (5) |Gray et aL] ( |2006^ . 



Table 2. Number of measurement, time span, and weighted 
r.m.s. of the residuals around the one planet solutions for 
HD 63765 for coRALiE-98, coralie-07, and harps. 



Instrument 



A',,, 



Span (T(o-c) 
[years] [ms"'] 



CORALIE-98 38 
CORALIE-07 18 
HARPS 52 



7.0 9.6 

3.18 4.72 

6.3 3.27 



A strong peak is present in the periodogram at ^56 day (see 
Fig. [3]l, which is identified as a Neptune-like planet in a circu- 
lar orbit by the genetic algorithm. We investigated whether the 
presence of stellar spots could mimic the short amplitude P - 
55.81 day Keplerian motion and explain such a laige x^- - 2.46 
since the derived orbital period - from activity index - is 34.7 
days. The bissector inverse slope, as well as the log [RfffA (see 
Fig.[T]i, do not present any variation at the period of P=55.81 day, 
while the measured radial velocities show a clear Keplerian sig- 
nature. The detection of the 55.81-day-period Neptunian planet 
is therefore confirmed. 

The residuals' dispersion is fairly large (cr - 4.6 m.s') and 
shows a clear low frequency oscillation at Pa!500 day that mim- 
ics the signature of a second Neptunian planet. However, this 
signal is correlated with the log [RfjiA activity index (see Fig. 
HI), ruling out the detection of a second long period planet in the 
system. Resulting orbital parameters are listed in Table |3] 



3.4. A 1 0-day-period Netpunian planet around HD 125595 

We have been observing HD 125595 with HARPS since May 
2004. Altogether, we have obtained 1 17 radial-velocity measure- 
ments with a typical signal-to-noise ratio of ^ 78 (per pixel at 
550 nm) and a mean measurement uncertainty (including pho- 
ton noise and calibration errors) of 1 ms"'. 

A clear « 10-day periodic signal can be seen in the peri- 
odogram (v=0. 1 day"') with its replica on both sides of the one- 
day period alias (see Fig. l5]l. A single-planet Keplerian model 
was therefore adjusted to the data leading to a P = 9.67 days 
circular orbit of semi-amplitude K - 4.79 ms"' (see Fig. pb, 
which corresponds to a planetary companion of minimum mass 
m sin (/) - 0.77 MNept- Corresponding orbital elements are listed 
in Table [3] The ~ 3 ms"' dispersion in the residuals is explained 
by a significant amount of energy in the periodogram of the 
residuals around the orbital solution at a period of 37 days (see 
Fig.|6]). This P ~ 37-day energy peak is also present in the peri- 
odograms of the CCF FHWM (resp. log(/?^^j and bisector in- 
verse slope and is close enough to the rotation period of the star 
to asses that the observed dispersion in the radial velocities is 
probably induced by activity. 



3.5. A long-period massive planet around HIP 70849 

HIP 70849 has been observed with HARPS since April 2006. 
Altogether, we gathered 18 radial-velocity measurements with a 
typical signal-to-noise ratio of » 49 (per pixel at 550 nm) and 
a mean measurement uncertainty (including photon noise and 
calibration errors) of 2.2 ms"'. A low frequency radial velocity 
variation with a large amplitude (x 200ms"') is clearly seen in 
the data (see Fig. iTJi. Even though the orbital phase is not cor- 
rectly covered by our observations, we were lucky enough to 
obtain radial velocity measurements close to the date of passage 
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Fig. 3. Top: Generalized Lomb-Scargle periodogram of the radial-velocity measurements for both HD 63765 (left) and HD 104067 
(right). Middle:Radial-velocity measurements as a function of Julian Date obtained with CORALIE (blue and orange dots) and 
HARPS (red dots). HARPS' only phase-folded radial velocities are also displayed (bottom). The best Keplerian, one-planet-solution 
is displayed as a dark curve whose corresponding orbital elements are listed in Table [2] 



through periastron. Running our genetic algorithm followed by 
MCMC simulations with the simplest model (1 planet, no radial 
velocity drift) allowed us to bring some constraints on the orbit 
as shown in Fig. [8] Indeed, at the end of the simulation the chains 
remained within bounding values with a period between 5 and 90 
years, an eccentricity between 0.4 and 0.98 and a minimum mass 
between 3.5 and 15 Mjup. It should be noted that solutions with 
much longer periods and very high eccentricities exist but could 
not be probed by our MCMC simulations in a reasonable amount 
of time. However, the corresponding minimum mass of such so- 
lutions remains within the substellar domain. The addition of a 
linear drift to the model decreases the planet eccentricity, as well 
as its radial velocity semi-amplitude, which as a result, decreases 
the mass of the detected object. We also looked at the Hipparcos 
data to check whether the astrometric signature of the compan- 
ion (a > 500juarcsec) is detected, to determine the real nature 
of this object as for HD 5388 b ( Sahl mann et al.|201 la i. In prin- 
ciple, astrometric motions larger that 1 mas could be detected 



in the Hipparcos intermediate data as shown by Sahlmann et al 



( 2011b| l, provided that the RV orbital elements are determined 
well and that the period is close to the duration of the Hipparcos 
mission. In that particular case, the long period of the planet 
combined with the poor determination of the orbital elements 
prevents to disentangle the companion signature from Hipparcos 
proper motions. However, future space missions such as GAIA 
( [Penyman et al.|2001| l or JMAPS (Hennessy et al.|2010[ ) will be 
able to characterize the nature of this object better. 



4. Conclusion 

We have reported the discovery of four extrasolar planet 
candidates orbiting moderately active stars - HD 63765, 
HIP 70849, HD 104067, and HD 125595 - discovered with 
the HARPS Echelle spectrograph mounted on the 3.6-m ESO 
telescope located at La Silla Observatory. Two of those planets, 
HD 125595 b and HD 104067 b, are Neptune-mass-planets 
while the other two - HD 63765, HIP 70849 - are clearly gas 
giants. 
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Table 3. Best Keplerian orbital solutions for HD 63765,HD 104067, HD 105595, and HIP 70849. 



Parameters 




HD 63765 b 


HD 104067 b 


HD 125595 b 


HIP70849 b 


r 


[kms-'] 


22.5843 ± 0.0024 


15.1243 ±0.0009 


5.2637±0.00054 


0.099±0.032 


AV(C07 - C98) 


[ms-'] 


-3.9 ± 2.4 








AV{H - C98) 


[ms-'] 


39.4±2.0 








'.vg 


[BJD-2.45 10«] 


4224.7731 


4001.5984 


4439.9245 


4564.01764 


p 


[days] 


358.0 ±1.0 


55.806 ± 0.049 


9.6737 ± 0.0039 


5-90 [years] 


K 


[ms-'] 


20.9 ±1.3 


11.56 ±0.75 


4.79 ± 0.47 


70-400 


e 




0.240 ± 0.043 


0.00 


0.00 


0.47-0.96 


O) 


[deg] 


122 ± 13 


0.00 


0.00 


- 


To 


[BJD-2.45 10'^] 


4404 ±11 


4043.15 ±0.56 


4435.23±0.15 


- 


-iCfavg) 


[deg] 


302.1 ±3.6 


91.9 ±3.6 


175±5 


- 


A(favg) 


[deg] 


302.3 ± 7.0 


91.9 ±3.6 


175±5 


- 


ai smi 


[10-^AU] 


667.8 ± 40 


59 ± 3.7 


4.2 ± 0.41 


11.7-661.10^ 


f(m) 


[lO-'^Mo] 


313 ±57 


9.0 ± 1.7 


0.111 ±0.031 


95 - 428 . 10^ 


m s'mi 




0.64 ± 0.05 Mj„p 


3.44 ± 0.25 MNep, 


0.766 ±0.078 MNep, 


3-15 Mj„p 


a 


[AU] 


0.940 ±0.016 


0.2643 ± 0.0045 


0.0809 ± 0.0014 


4.5-36 


N 




108 


88 


117 


21 


Time span 


[years] 


10.77 


6.22 


5.68 


4.22 


o-(o-c) 


[ms-i] 


3.41 


4.60 


3.22 


3.94 - 3.68 




X7 




1.30 ±0.16 


2.46 ± 0.24 


1.61 ±0.17 


1.06 ±0.38 



Notes. 4vg is the mean date of the observations, /l(?avg) is the mean longitude at t 
dispersion of the residuals. 



g, A(4vg) is the true longitude at t 
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Fig. 4. Correlation bewtween the activity index and the radial 
velocity residuals as function of time for HD 104067. Pearson's 
r = 0.59. 



- //D 725595 Z?: According to Mordasini et al. (20091, the 



planet is too massive {nip - Q.ll MNept) to have formed "in situ" 
since there are not enough planetesimals at 0.08 AU with typical 
disk profiles. Because its parent star has solar metallicity, the 



planet probably formed beyond the iceline and then migrated 

toward the star through tidal interaction with the gaseous 

disk. The planet is not massive enough to trigger runaway gas 

accretion. 

- HD 1 04067 b .The existence of such a planet, however, is 

hardly explained by planetary formation and evolution scenarios 



(Ida & Lin 2005 



Mordasini et al. 2009 1. Indeed, the planet 



IS massive enough (m/,=3.4]V[Nept) to have triggered runaway 
gas accretion before migrating toward the star. How could the 
detection of such a planet be explained? Are we biased by the 
observing techniques since HD 104067 b could be a gas giant 
provided that its inclination is lower that 30 deg? This question 
will be answered when such an object is detected in transiting 
surveys or when a statistically reliable sample of similar object 
is detected by radial velocity surveys. 

- HD 63765 fe:This is a one-year-period Saturn-mass planet with 
a significant eccentricity. There is nothing peculiar about this 
object but its detection is useful to build a statistically robust 
sample of giant planets properties (orbital elements and mass 
distributions, host star properties). 

- ///f 70iS49 fo:Interesting is the presence of a massive giant 
planet around the 0.63 solar-mass-star HIP 70849. Even though 
characterized as a K7 dwarf, HIP 70849 is close enough to 
the very low-mass star domain to present similar planetary 
formation scenarios . According to [La ughlin et al. ( 2004|, [Idal 
|& Lin| P005| l, and [Mordasini et al.| ( |2009| ), the formation of 
Jupiter-mass planets orbiting very low-mass stars is seriously 
inhibited at all separations. The presence of this 3-15 Mjup 
planet in orbit to HIP 70849 brings to four the number of 
companions more massive than 2Mjup in orbit around a very 
low mass star and therefore questions the validity of planets 
formation scenarios around very low -mass stars. Possible hints 
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Fig. 5. The top figure is the periodogram of HARPS radial veloc- 
ity measurements of HD125595. A clear signature is seen at fre- 
quency v=0.9 day ' with its replica on both sides of the one-day 
alias. The bottom figure represents the phase-folded radial veloc- 
ity measurements with the best Keplerian, one-planet-solution 
represented as a black curve. 

of the existence of such objects could be that the mass of their 
protoplanetary disk did not scale with the mass of the star or, 
otherwise, that the planet started to form much earUer than 
usually thought. 
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Fig. 7. Radial-velocity measurements as a function of Julian 
Date obtained with HARPS (red dots) on HIP 70849. One of the 
orbital solution is overplotted as a black line. 
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Fig. 8. The top diagram represents the joint probability density 
function of the orbital period and the eccentricity. The red (resp. 
blue, cyan) line corresponds to the l(resp. 2, 3)cr confidence 
interval. The bottom diagram represents the minimum mass dis- 
tribution of the planetary companions with the same confidence 
intervals. 99.73% of the solutions have a period between 5 and 
90 years, an eccentricity between 0.47 and 0.95 and a minimum 
mass between 3.5 and 15 Mjup. 
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